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Abstract 

This paper aims at designing an on-board beam generation process for multibeam satellite systems 
with the goal of reducing the traffic at the feeder link. Full frequency reuse among beams is considered 
and the beamforming at the satellite is designed for supporting interference mitigation techniques. In 
addition, in order to reduce the payload cost and complexity, this on-board processing is assumed to be 
constant and the same for forward and return link transmissions. To meet all these requirements a novel 
robust minimum mean square error (MMSE) optimization is conceived. The benefits of the considered 
scheme are evaluated with respect to the current approaches both analytically and numerically. Indeed, 
we show that with the DVB-RCS and DVB-S2 standards, our proposal allows to increase the total 
throughput within a range between 6% and 15% with respect to other on-board processing techniques 
in the return and forward link, respectively. Furthermore, the proposed solution presents an implicit 
feeder link bandwidth reduction with respect to the on-ground beam generation process. 
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I. Introduction 

The increasing demand for fixed broadband data services is an opportunity for satellite indus¬ 
tries to target new markets apart from the well-known current ones (i.e. broadcast broadband, 
emergency communications, ...)• In order to cope with higher data traffic demands, satellite 
system designers are looking for advanced satellite communication architectures. In this context, 
the use of multiple beams has recently received a lot of attention as a key enabler of next 
generation high throughput satellite systems. These systems rely on implementing a large number 
of beams instead of a single (global) beam in the coverage area. This is beneficial since each 
beam can have a larger antenna gain-to-noise temperature than in the single beam case and the 
available spectrum can be reused among spatially separated beams. Furthermore, whenever the 
satellite systems delivers broadband interactive traffic, the multibeam architecture can support 
different modulations and code rates to each beam depending on the user link quality, leading 
to a high increase of the overall system throughput. 

The implementation of multibeam satellite systems is currently being investigated in order to 
increase the overall spectral efficiency while keeping the payload complexity affordable. One 
of the main challenges is how to deal with the large spectral demands of the feeder link (i.e. 
the bidirectional communication link between the satellite and the service provider), whose 
bandwidth requirements increase exponentially as it aggregates the traffic of all users. Recently, 
some techniques have appeared in order to reduce the feeder link spectrum requirements. There 
is a current tendency for moving the feeder from the Ka band to the Q/V band, where there are 
larger available bandwidths [[I]|. Unfortunately, in these frequencies the fading is extremely large 
and more advanced transmitting diversity techniques are needed. 

Another option is the use of multiple gateways, which might be adequate in order to reduce the 
feeder link spectral requirements as they can be equipped with very directive antennas and exploit 
the spatial diversity while sharing all available spectrum (H, [l3]|. Nevertheless, the deployment 
of several gateways increases the cost of the system and; moreover, the interference mitigation 
techniques suffer from certain degradation lIH, [15]|. This is due to the fact that the processing 
must be separated in isolated processing units. 

In contrast to the aforementioned feeder link traffic reduction techniques, this paper focuses 
on the on-board beam generation process. This promising solution keeps certain processing in 
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the payload so that the amount of required signals from the feeder link are severally reduced. 
In this way, the satellite does not act in transparent mode and it carries out some processing, 
leading to a high reduction of the feeder link bandwidth requirements. Specifically, while the 
on-ground beamforming requires a feeder link bandwidth of 




feeder link on-ground 


= 


( 1 ) 


where N is the number of feed elementqj and i?user is the user bandwidth; the on-board beam¬ 
forming only requires 


B 


feeder link on-board 




( 2 ) 


where K is the number of users. For this work we will consider multiple-feed-per-beam archi¬ 
tecture where N > K. Note that, in contrast to single-feed-per-beam architectures {N = K), in 
multiple-feed-per-beam architectures beamforming scan losses are negligible ||6l. A more detailed 
description of the beam process is presented in [|71|. 

Apart from the feeder link challenge, multibeam satellite systems require a large capacity in the 
access network. As a matter of fact, in the generated radiation pattern on Earth, adjacent beams 
create high levels of interference and, therefore, a carefully planned power and frequency reuse 
among beams must be employed to cope with this increased level of interference. Consequently, 
beams with adjacent footprint currently operate in different frequency bands or polarizations. 
In this context, the number of colors Nc is the essential parameter, which corresponds to the 
number of disjoint frequency bands and polarizations employed on the coverage area > 1). 
In fact, the lower the number of the higher the overall system bandwidth will be and the 
higher the interference power levels will be generated. 

In order to increase the available bandwidth yet maintaining a low multiuser interference, a 
promising technique is to use full frequency reuse pattern {Nf. = 1) and resort to interference 
mitigation techniques. In this way, signals can be precoded and detected before being trans¬ 
mitted and received in order to reduce inter-beam interference lO. As a result, a considerable 
improvement of the achievable spectral efficiency can be obtained. To this end, more advanced 
interference mitigation techniques as precoding in the forward link and multiuser detection or 


*The input signals of the antenna array feed assembly located in the payload. 
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filtering in the return link have been considered in past studies of the European Space Agency 

(ESA) m, a. 

Since interference mitigation techniques require large computational resources, they must be 
carried out on ground. Indeed, larger efficiencies are obtained if not only the precoding and 
detection are done on ground, but also the beam generation process, as more flexible processing 
units are available. In other words, if the beamforming is kept fixed on the payload, there is 
a performance loss compared to the spectral efficiencies obtained by on ground beamforming 
m, EH- However, if the satellite does not perform any beam processing, the feeder link needs 
a large amount of spectral resources in order to transmit all the user signals and beamforming 
weights. Consequently, even though certain degradation is expected with respect to the on-ground 
operation (i.e. beam generation, precoding and detection are done in the terrestrial segment), in 
the present work we propose to optimize the on-board beam generation process so that the 
achievable rates do not severally decrease due to the on-board beam generation and the feeder 
link traffic is kept low. 

Concretely, this paper focuses on obtaining an optimal on-board beam generation when linear 
minimum mean square error (EMMSE) precoding technique in the forward link and EMMSE 
detection procedure in the return link are used as interference mitigation techniques. This study 
foresees the presence of a non-channel-adaptive (fixed) on-board beam processing scheme in 
order to keep payload complexity low. Thus, the problem becomes more difficult in the presence 
of this fixed process in the payload. In order to deal with this problem, we use a robust 
optimization framework so that a fixed beam generation can be obtained despite user link channel 
variation. 

Eurthermore, the design for both the forward and return links results the same, which makes it 
appropriate for the future multibeam satellite systems since it is expected that the same reflector 
is employed at the return and forward links. Note that the variability of the channel is mainly 
due to the change of position of the users in consecutive time instants. Numerical simulations 
show the benefit of our method, which in some scenarios can increase the spectral efficiency 
over the 6% and 15% for return and forward links, respectively, if the DVB-S2 and DVB-RCS 
modulation and coding parameters (modcods) are used. 

To the best of the authors knowledge, this is the first time the problem of on-board beam 
generation process is treated not only in the forward but also in the return link. In contrast to 
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our preliminar work IfT^ . where only the forward link was examined, in this paper we focus our 
attention to the joint forward and return link optimization. In addition, a novel and better robust 
design is presented based on a tighter upper bound of the optimization problem. This new scheme 
is conceived considering a first order perturbation approach. Finally, several detailed evaluations 
are presented that validate our contribution in detail. 

The rest of the paper is organized as follows: Section II presents the signal model. A brief 
introduction of the beam generation process and the problem characteristics are described in 
section III. Section IV presents the novel solution that the paper proposes. Section V presents a 
novel robust scheme based on a first order perturbation analysis. Section VI contains a summary 
of the simulation results, and eventually the conclusions are given in section VI. 

Notation: Throughout this paper, the following notations are adopted. Boldface upper-case 
letters denote matrices and boldface lower-case letters refer to column vectors. (.)^, (.)^, (.)* 
and (.)+ denote a Hermitian transpose, transpose, conjugate and diagonal (with positive diagonal 
elements ) matrix, respectively. Itv builds N x N identity matrix and Okxn refers to an all-zero 
matrix of size KxN. If A is aNxN matrix, Ai;^ refers to taking the K first rows of the 
matrix A. (X)jj represents the (f-th, j-th) element of matrix X. If B is a V x matrix, A < B 
implies A — B is semidefinite negative, a -< b means vector a majorizes vector b. Finally, E{.} 
and ||.|| refer to the expected value operator and the Frobenius norm, respectively.- denotes the 
matrix Hadamart product. 


II. Signal Model 

Let us consider a multibeam satellite communication system, where a single geosynchronous 
satellite with multibeam coverage provides fixed broadband services to a large set of users. To 
this end, the satellite is equipped with an array fed reflector antenna whose number of feeds is 
denoted by N. The coverage area is divided into K beams, with 

K <N, (3) 

and the users are assumed to be uniformly distributed within the beams. By employing a time 
division multiplexing access (TDMA) scheme, at each time instant the gateway is serving a total 
of K single antenna users (i.e. exactly one user per beam), and it is transmitting (receiving) 
information to (from) the same number of the users through the satellite in the forward (return) 
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link. Note that in return link satellite eommunieations generally operate in a multi-frequeney 
TDMA (MF-TDMA) so that different users of the same beam might be alloeated to different 
sub-bands. For the sake of simplieity and without loss of generality, the rest of the paper eonsiders 
TDMA for the return link. Remarkably, the eoneeived teehnique ean be aeeommodated to the 
multi-band eommunieation by replieating the linear proeessing at eaeh band due to the frequeney 
flatness of the ehannel response. 

The satellite is assumed to linearly eonvert a set of N on-board feed signals into the K 
feeder link signals whieh are transmitted to the gateway in a frequeney multiplexed fashion. 
Reeiproeally, in the forward link, the same linear proeessing strategy is used to eonstruet the N 
feed signals from the K feeder link signals. 

Moreover, sinee a high throughput system is targeted, full frequeney reuse among beams 
is assumed so that all beams ean share the same frequeney resourees. The user link is the 
eommunieation bottleneek of the whole system. The feeder link is assumed perfeetly ealibrated 
and noiseless. Figure [T| summarizes the transmission bloek diagram. 

In the following, the mathematieal expressions of the signal model in both the return and 
forward links are deseribed. 

A. Return Link 

As stated above, K denotes the number of users and N is the number of on-board feeds. 
Then, the eorresponding reeeived signal at the gateway ean be modelled as 

Yru = ^BHs Bn, (4) 

where ym = [vrl,!, ■■■■,yRL,KY is a K x 1 veetor eontaining the staek of reeeived signals at 
the gateway. The K x 1 veetor s is the staek of the transmitted independent signals by all users 
sueh that E{ss^} = I^. Note that, throughout the paper the subseript rl is used to refer the 
return link while fl will denote the forward link. The eonstant (3 denotes Equivalent Isotropie 
Radiated Power (EIRP), whieh is referred to the user terminal transmit power and we assume 
to be the same for all the users. 

In order to radiate the multibeam pattern, the satellite payload is equipped with a smart 
antenna system (generally an array fed refleetor) eoined as on-board beam generation proeess. 
This system eonstruets the beam pattern for transmitting and reeeiving data from the eoverage 
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Fig. 1. Multibeam satellite system with on-board beam generation process. The precoding and detection procedures are done 
on ground. On the contrary, the beam generation process is carried out at the payload and it is assumed to be constant and the 
same for forward and return links. 


area. Mathematically, the effect of this beamforming appears as the rectangular KxN fai matrix 

B. 

The N X 1 vector n accounts for the zero mean Additive White Gaussian Noise (AWGN). 
We assume unit variance Gaussian noise samples such that 

E{Bnn^B^} = BB^. (5) 

For radio-frequency design convenience, we will assume that B is orthonormal so that the feed 
signals are decoupled at the payload (BB^ = \k) ■ Matrix H is the overall N x K user link 
channel matrix whose element hjj presents the aggregate gain of the link between the z-th satellite 
feed and the j-th user (in the j-th beam). This channel can be decomposed as follows: 

H = GD, (6) 


where: 
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• G is a X K matrix that models the feed radiation pattern, the on-board attenuation and 
path losses. It is responsible for the interferenee among users. We assume the elements of 
G are normalized so that they have unit varianee. 

. D is assumed to be a if x if diagonal matrix which takes into account the atmospheric 
fading in the user link. 

Let us describe more exhaustively matrix G. Its {k^n)-\h entry can be described as follows 


(G) 


k^n 


Cl R^kn 


(7) 


with dk the distance between the /c-th user terminal and the satellite. A is the carrier wavelength, 
it's is the Boltzmann constant, Bw is the carrier bandwidth, G\ the user terminal receive antenna 
gain, and Tji the receiver noise temperature. The term akn refers to the gain from the n-th feed 
to the k-\h user. It is important to mention that the G matrix has been normalized to the receiver 
noise term. The reader can refer to [fTTlI for a more detailed description of the channel model. 


B. Forward Link 

Analogously to the return link, the signal model of the forward link becomes 

Yfl = + w, (8) 

where K x 1 vector ypu is the stack of received signals at each user terminal, and x is a 
K X 1 vector that contains the stack of transmitted symbols. Remarkably, in general wireless 
communication systems, the channel reciprocity does not hold as uplink and downlink operate 
in disjoint frequency bands. However, considering our channel modelling , the channel matrix in 
the forward link differs from the return link in the path loss, feed gain and atmospheric fading. 
As a result, a scaling factor 7 can model the different frequency operations. As explained later 
on, this rescaling factor does not influence the proposed optimization and; therefore, it can be 
set 7 = 1 . 

Similarly as in the return link, w is a iT x 1 vector that represents the independent and 
identically distributed zero mean Gaussian random noise with unit variance such that 


£'{ww'^} = 


( 9 ) 
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Evidently, B does not influence in the forward link noise covariance matrix. We assume the 
following average available power constraint: 

trace (xx^) < Pfl, (10) 

where Ppu denotes the total transmit power in the forward link. Note that the transmit power 
constraint is set without considering the beam generation process B. This is because the power 
allocation mechanism is located before the array fed reflector system. In addition, it is assumed 
that the feeds can share the available transmit power. This can be implemented with flexible 
power amplifiers as described in lfT3l . 

Now, we proceed to jointly optimize matrix B so that the overall system performance is 
improved. It is important to remark that B must be the same for both the optimization of the 
return and forward links in order to reduce the payload cost. In addition, this matrix needs to 
be constant in order to keep the payload complexity low and minimize the feeder link spectral 
resources. 


III. Problem Formulation 

Let us assume that the gateway has perfect Channel State Information (CSI) and uses LMMSE 
as described in m for precoding in the forward link and LMMSE filtering for multiuser 
detection in the return link. These techniques have been pointed out as efficient methods due 
to both its interference rejection capabilities and fairness among beams while preserving a low 
computational complexity ifTSl . 

This work resorts to the minimization of the trace of the MSB matrix both at the forward and 
return links that results from the use of LMMSE precoding and detection. Let us briefly outline 
the overall mathematical derivation: 

1) First, the MSB matrix of the return link is computed assuming LMMSE detection. 

2) Second, the MSB matrix of the forward link is computed assuming LMMSE precoding. 

3) Third, an upper bound of the MSB minimization in the return link is presented. 

4) Finally, a novel robust beam generation process in the return link, which considers the 
aforementioned upper bound is obtained. For the forward link case, the optimal design 
yields to the same solution as it is described. 
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Remarkably, the design of the optimal B is imposed to be non ehannel dependent. We show 
that the optimal B in the forward and return links results to be the same; thus, fulfilling one of 
the eonstraints of the system. 

A. Return Link 

As a first step, let us define as the LMMSE filter that deteets K reeeived signals at the 
gateway such that s = composed by s* which denotes the i-th element of the detected 

signal (for Athe user) in the gateway. In this context, the MSE of Ath user is achieved as follows 

MSERL,, = E{|si-Si|2}, (11) 

where s* represents the Ath element of transmit signal vector (for Athe user) for a total of K 
users such that s = (si,..., skY . 

It is well known that the mathematical expression of LMMSE filter becomes 

= {\k + /3H^B^BH) H^B^, (12) 

and the MSE matrix after the use of this filter is 

MSErl = (^Ix + /3H^B^ (BB^)"^ Bh) ^. (13) 

Without loss of generality, we restrict B to be orthonormal such that BB^ = \k- The sum 
of MSE in the return link is defined as 

SMSErl = trace (l^ + /5H^B^BH) . (14) 

Now, let us assume for a moment that B can be channel adaptive (i.e the payload can modify 
B depending on the channel variations) . Then, the corresponding problem is formulated as 

(Ix + /9H^B^BH)~') (15) 

s.t. BB^ = \k- 

It is important to remark that the authors in ifTOll showed that the presence of B increases the 
SMSErl in the gateway. Mathematically, 

trace (^{Ik + /3H^B^BH)"^) > trace . (16) 


min trace 

B 
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Indeed, in IfTOll it was shown that with the following Singular Value Deeomposition (SVD) of 
the ehannel H = U$V^, an optimal design of B ean be worked out as 

B = (17) 

where denotes the K first rows of the matrix U^. In faet, it ean be easily seen that this 
partieular solution reaehes equality in (fT^ and; thus, minimizes the SMSErl- 

In the present work, B is assumed to be non-ehannel adaptive, therefore, the design of B in 
(fTTl) eannot be eonsidered. Even though the ehannel appears to be variable at eaeh realization, 
we aim at finding the best possible non-ehannel adaptive design of B. In this eontext, let us 
deeompose the ehannel as follows 

H = H + A, (18) 


where: 

• H represents the mean value of the ehannel. 

• A models the differenee between the aetual value of the ehannel and its mean. It indieates 
the variability of the ehannel in eonseeutive time instants as already explained in seetion I. 

We assume that the aetual ehannel H lies in the neighbourhood of a nominal ehannel H that 
is known to the gateway. In partieular, we eonsider that H belongs to the uneertainty region 
7( = {H: ||H — H|| < a} whieh is an sphere eentered at H with the radius a. 

Interestingly, the ehannel model in (fTSl) resembles the modeling of a MIMO system with 
imperfeet CSI at the transmitter whieh has been solved as a worst ease optimization problem 
in ifT^ - lfT^ . With this perspeetive for the return link, the worst ease robust design is proposed, 
whieh leads to a maximin or minimax formulation: 


mm 

B 


max 

A 


traee + 


-1 


(19) 


s.t. 


BB^ = Ik. 


Prior to obtaining the solution of (fT9l) . let us foeus on the forward link optimization problem, 
whieh is similarly derived. 


B. Forward Link 

In the forward link, the zero foreing preeoding with a regularized inversion is assumed [[T4l . 
In this ease, the linear preeoding is expressed as 


X = Tc, 


(20) 
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where T is the K x K preeoding matrix at the gateway and c is the K x 1 transmit symbol 
veetor at all feeds sueh that E{cc^} = Ik- In this eontext, the eorresponding preeoding matrix 
T is expressed as 


T = ^B*H* ( -^Ik 

TFL 




( 21 ) 


where the value of the eonstant p has to eomply with the forward link power eonstraint as 
follows 

traee (TT^) < PpL- (22) 

This partieular kind of preeoder is used to find an optimal balanee between aehieving signal 
gain and limiting the multiuser interferenee. Similar to the return link, MSEpL.i is defined as 

MSE,.i,i = E{|ci-Ci|2}, (23) 

where MSEpu,! refers to the MSE received by Pth user. Similarly, c = (ci, ...,ck)^ and c = 
(y^)“^yFL = (ci, ...,ck)'^ are the transmitted and received signals for K users, respectively. In 
this context, c* represents the transmitted signal for Pthe user and c* denotes the signal received 
by user Pth. The MSE matrix in the forward link can be calculated as follows 


MSEpT = 


which can be rewritten as in 


MSEpL = ^ ( fH^B^B*B^B*H* + -^Ik 
Pvh \ V -Pel 


e{((Vp) Vfl-c) {{,/p) Vfl-c)^|, 


H^B^B*H* + -/-Ik 
Pel 


K 


-2' 


As in the return link, we concentrate our efforts to minimize the sum of MSE, this is 

SMSEpL = trace(MSEpL), 


(24) 


(25) 


(26) 


where, recalling that BB^ = Ik and we consider the following property, trace(A) = trace(A 


T\ 


where A is a square matrix, then we have that 

K 


SMSEfl = -/-trace I ( H^B^BH + -^Ik 
Pel \ V Pel 


K 


- 1 ' 


(27) 


The worst case optimization problem thanks to the channel decomposition in (fTSl) can be 
formulated as follows 

(28) 


min max trace I ( H^B^BH + ^^1k 
BA \ V PpL 


K 
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s.t. BB^ = Ik. 

Note that the return link optimization (fT9l) and the forward link one (l28l) are the same exeept 
for a scalar value. In next section we show that both lead to the same optimal design; thus 
confirming a natural uplink downlink physical duality. 

IV. B OPTIMIZATION 

This section tackles with the main objective of this paper. An optimally designed B for 
problems (fT^ and (1281) is presented. Two main steps are followed. The first step provides a brief 
description of an upper bound for the SMSE. The second step proposes a design for B such 
that it minimizes the proposed SMSE upper-bound obtained in the first step. The design is done 
for the return link and extended to the forward link. 

Prior to presenting the optimal design, we need to introduce the next lemma. 

Lemma 1. Assuming an arbitrary square matrix A, the next equation holds 

trace ^(Ijt + AA^) = trace + A^A) . (29) 

Proof. It is a direct consequence of inversion matrix lemma. □ 

By considering A = v(5BH, the SMSErl in problem (fT9l) can be rewritten as 

trace + /5BZB^)"^) , (30) 

where Z = = HH^ + HA^ + AH^ + A A^ is a V x matrix. We propose an upper 

bound of SMSErl as follows 

Theorem 1: The SMSErl is upper bounded by 

trace (^{Ik + < trace {^K + (31) 

where 

Z ^ U(S - enlNyV^, (32) 

so that HH^ = USU^ is the eigen-decomposition of matrix HH^. The scalar value ch is 
defined as 

eu = (33) 

where 5max{C>) denotes the maximum singular value of C matrix. 
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Proof. See [ETI see.7.3.1]. □ 

As a result, a worst-ease SMSErl ean be obtained in practice by using the lower bound Z in 
lieu of Z. However, it is important to mention that some values of a lead to infeasible MSErl 
solutions, that is, for a large value of a the matrix (l32l) might become low rank since O'*" operator 
delivers 0 whenever the diagonal entry is nonpositive. In order to avoid this circumstance, the 
value of a has to be checked and, if necessary, decreased so that the feasibility condition of the 
problem (fT9l) is hold. 

In order to obtain a robust design, the target is to minimize the proposed upper-bound of 
SMSEjil in (I3TI) instead of (fT9l) . In this case, the corresponding problem is formulated as 


min 

B 


trace 


'-K 


/3BZB 


H 



s.t. BB^ = Ik. 


(34) 


The solution to this optimization problem is described in the next theorem. 

Theorem 2: Let B and be two matrices of size K x N and N x N, respectively. Then, 
the upper bound of SMSE is minimized if B is selected as the first K rows of the matrix L^, 
that is 

B* = (35) 

where B* denotes the optimal design ofB. 

Proof. See Appendix A. □ 

Remark: It is important to mention that the derivation of theorem 2 differs to theorem I 
in |fT9l . The main difference relays on the constraint since in lfT9l a total power constraint is 
considered 

trace (BB^) < P, (36) 

where as this paper assumes 

BB^ = liT, (37) 


which involves further mathematical developments as described in Appendix A. 
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Before starting with the forward link case, let us remark that B* only needs statistical channel 
knowledge in order to be computed. Moreover, its design does not depend on ot. Indeed, the 
value of a affects only on the resulting SMS Erl- This is due to the optimization of an upper 
bound of the problem instead of the problem itself. Now, let us proceed with the forward link 
optimization. 

In the forward link the optimization problem can be formulated as follows 


mm trace I I BZB^ + 

B \ V Pfl 


(38) 


s.t. 


BB 


u 


= lx- 


In can be observed that the optimal solution of (l3^ is (1351) . The sketch of the proof is similar 
to the one presented previously for the return link and; thus, we only comment it. The idea 
is to check whether the term does not influence the optimal value of (l38l) which can be 
easily observed in appendix A. Consequently, the scaling factor due to the channel variations 7 
does not influence the optimization, either. Remarkably, this derivation is different from the one 
presented in our preliminary work in ifT^ . because this paper so encompasses the forward and 
return link optimizations. 

Note that the robust beamforming design has the same eigenvectors as the nominal channel 
matrix HH^. In other words, the presented robust design only considers eigenvalue variations 
due to the different user positions. In the next section, the impact on the eigenvectors is analysed. 

V. First Order Perturbation Analysis 

As discussed in the previous sections, the underlying optimization problem (|2^ shall be lower 
bounded in order to obtain a closed-form solution. This is done by means of considering upper 
bounds of Z. 

Indeed, the proposed perturbation model can be described as 

Z = (tj, + AU,) (S, + AS,) (U, + AU,)"" + (U„ + AU„) (S„ + AS„) (U„ + AU,)"", 

(39) 

where the U denotes the matrix containing the eigenvectors and S is a diagonal matrix which 
contains the eigenvalues. Subindex s denotes the non-zero signal space whereas n the signal 


SUBMITTED 


16 


space that is spanned by the zero valued eigenvalues (i.e. the null spaee of Z). All AUs, 
ASs, AU„, AS„ are generated by a perturbed version of Z: 

Z = Z + AZ, (40) 

where 

Z = HH^, (41) 

and 

AZ = HA^ + AH^ + AA^. (42) 

Under this eontext, U denotes the eigenveetor of the nominal matrix Z whereas S a matrix 
eontaining its eigenvalues. The other matriees with the A- prefix denote the eorresponding 
perturbation matriees. 

The previous seetion has implieitly eonsidered two assumptions. First, it has been assumed 
that the ehannel variations do not modify the dimension of the null spaee so that AS„ remains 
as a zero matrix. Seeond, it has been assumed that AU^, = 0 ,whieh might not be true in eertain 
eases lf20ll . The aim of this seetion is to eonsider the effeet of this later perturbation in order to 
obtain a tighter upper bound of Z than the presented in the previous seetion. Remarkably, the 
following inequality holds 

Z > Z > Z, (43) 

where Z only eonsiders perturbations at the eigenvalues whereas Z eonsiders both perturbations 
at both eigenvalues and eigenveetors (AU^). Next theorem provides an approximate solution 
whenever these both perturbations are eonsidered. 

Theorem 3 The beamforming matrix that optimizes the MSE upper bound when considering 
both eigenvector and eigenvalue perturbations is 

B'* = U(S-6 Hliv)^U^, (44) 

where 

U = tj, + (45) 


and 


R = D ■ (Uf U,S + SUf U,) , 


(46) 
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and the g, f-th entry o/D is 


A/-A,’ 


for f g and A/ for / = 1 ,..., denote the eigenvalues of 
Proof See Appendix B. 


(47) 


□ 


Note that for this ease, the eigenveetors of the beamforming matrix take a different value from 
the nominal matrix. In addition, the larger a the more different are the eigenveetors from the 
nominal channel matrix ones. 

As we have already seen, the beam generation process both on the forward and return links 
leads to the same matrix B, which is fixed. Now, it is time to compare the benefits of this design 
in front of the current beam processing deployments. 


VI. Simulation Results 

In order to show the performance of our proposal, this section presents a numerical evaluation 
of the conceived technique. Our baseline scenario is an array fed reflector antenna and matrix 
B that have been provided by ESA in the framework of a study on next generation multibeam 
satellite systems. The number of feeds is assumed to be = 155 and K = 100 beams that are 
covering the whole Europe area. 

Results have been averaged over a total of 1000 user link channel realizations. Note that, 
only atmospheric fading due to rain effect is considered in the user link channel and further 
refinements of the channel are neglected. This simple characterization is useful for the intended 
comparisons and it is a general practice in the evaluation of multibeam satellite systems. 

The randomness of the channel is due to the user positions which are assumed to be uniformly 
distributed within the beams. In addition, we will assume that each user employs all available 
spectrum and the atmospheric fading is modelled as in [f2T]| . 

Recall that, full frequency reuse among beams and noiseless feeder link have been considered 
in this work. In the sequel, we compute different performance metrics. Eirst, the SINR for each 
user after employing interference mitigation techniques among users is presented. Then, with that 
SINR value, the throughput is inferred according to DVB-RCS and DVB-S2 standards for the 
return and forward links, respectively ll22]| . If2^ . Eurthermore, the simulation results also provide 
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the associated Cumulative Distribution Function (CDF) of SINK which shows the availability of 
the user link. In this case, the instantaneous availability indicator for the /c-th user is given by 

= giSlNRk) (48) 


which is equal to 0 if the user link is unavailable (i.e, if the instantaneous SINK is lower than 
that required by the lowest modcod for the return link, i.e. SINR^ < 1.7dB, and for the forward 
link, i.e. SINR^ < —2.72dB ) and is equal to 1 otherwise. We also present the Shannon capacity^ 
obtained from the user SINR, 


C'shannon = log 2 (l + SINR), (49) 

and assuming that interference is treated as Gaussian noise. This measurement serves us to see 
the potential of our work independently of the satellite standard modulations and channel coding 
both for the forward and return links. 

Another performance metric to be considered is the fairness among beams. Note that this is 
of great interest for satellite operators where near to equal achievable data rates per beam are 
the target. For this purpose, we present the throughput index of dispersion, defined as 

(TtIi 

Index of Dispersion =-, (50) 

/^Th 

where ciTh and /ixh correspond to the variance and the mean of the user throughputs, respectively. 
This metric provides an indicator of how the data rates are dispersed with respect to the mean. 
The larger the index of dispersion is, the less the fairness the system achieves. 

For a best practice, as upper bound for the achievable rates we consider only on-ground 
processing at the gateway (i.e. no on-board processing) as it is described in ifTOll . From the 
return link point of view, the received signal dH), which is based on this on-ground scenario, is 
rewritten as 

Yrl = T^.ground (Hs + n), (51) 

where 

Wo„.ground = H (HH^ + (52) 

denotes the LMMSE detector filter at the gateway. Note that the linear processing is similar 
to (O but in this case it has been assumed that no beam processing is done. Considering the 


^Of course, we refer to the use of the Shannon formula instead of the channel capacity. 
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forward link, the received signal by the user terminals with this on-ground technique can be 
represented as 

ypL H Ton-ground^ (53) 

It is important to remark that although large data rates can be obtained if all the processing 
is carried out on ground, the required feeder link spectral resources severally increase, leading 
to an inefficient system. 

To sum up, in order to test the validity of the derived theoretical results in section IV, we 
compute the spectral efficiency of the following multibeam satellite system using precoding and 
detection algorithms for forward and return links respectively: 

• B based on a geographical reasoning (reference). 

• B* proposed by this study in (l35l) . 

• B* proposed by this study in (HD). 

. On ground processing (upper bound). 

In the sequel, the results are separated into two different subsections, return and forward link. 
In this context, the same fixed optimal design of on-board beamforming matrix is computed 
since this optimal design depends on the right eigen vector of channel average matrix, H. This 
is computed empirically considering the aforementioned 1000 channel user realizations. 

A. Return Link 

The return link operates at 30GHz, and is based on DVB-RCS standard lf22]| and we target a 
Packet Error Rate (PER) of 10“^. Eigure [2] depicts the evolution of the total average throughput 
(bits/symbol) as a function of the user EIRP (/3) for different scenarios. Although by means 
of using the DVB-RCS standard the obtained throughput gain is limited when the Shannon 
capacity is considered, higher gains are obtained with respect to the reference scenario . In other 
words, other modcods design would improve the benefits of the proposed technique with respect 
to the reference scenario. Note that the proposed robust design that consider the eigenvector 
perturbation improves the system throughput with respect to the design that only considers 
eigenvalue variations. Indeed, our proposal is approaching the upper bound of the on ground 
design. 

The corresponding availability probability is also provided in Eigure |3l In this case, our 
proposal also improves the reference scenario, leading to an increase of the system availability. 
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Fig. 2. Return link throughput values over different user EIRP (/3). 
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Fig. 3. Return link availability. 
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Fig. 4. Return link throughput index of dispersion. 



Fig. 5. Return link throughput with respect to channel variations. 


Remarkably, the fairness among beams is also improved as it is depieted in Figure 01 Lower 
values of dispersion index are obtained with our teehnique with respeet to the referenee design. 
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Finally, we study the impaet of the ehannel variations on the beam proeessing design. Bearing 
in mind that a in (l3^ determines this variation, we eompute this value and we present its 
eorresponding average throughput values in Figure [H The values of a are seleeted so that the 
feasibility of MSErl in (l3^ holds. It implies that 

(S — eHlAr)^ > 0 Vi = l,...,iV. (54) 

For a large value of a the matrix (l54l) might beeome semidefinite negative and; thus, ehanges the 
nature of the problem. In order to avoid this, a has to be eheeked so that the matrix (l54l) always 
remains semidefinite positive. It is observed that the larger a values, the less the throughput is 
obtained due to the ehannel mismateh. 

B. Forward Link 

The forward link is assumed to operate at 30GHz and is based on DVB-S2 standard with 
a PER of 10“®. Note that the working points were extrapolated from the PER eurves reported 
in the DVB-S2 guidelines doeument [|23]|. Based on lf23]l . it is possible to find a relationship 
between the required reeeived SINR and the speetral effieieney aehieved by DVB-S2 standard. 

The results are presented for the total bandwidth and as a funetion of the total available power 
denoted by Pel- Eigure [6] depiets the aehieved results of speetral effieieney and Eigure |7] shows 
the availability of the users in the forward link. Clearly, the proposed teehniques perform better 
than the benehmark system and again the robust design based on the eigenveetor perturbations 
behaves better than the one that only eonsiders the eigenvalues. 

The expeeted result of throughputs in Eigure [6] is justified by the availability in Eigure |7lln 
other words, the system with new proposed design of B is eloser to upper bound seenario than 
the referenee. Moreover, the impaet of ehannel variations ean be observed in Eigure HI It is 
elear that our proposal results in higher throughputs even when the ehannel variations are high. 
Remarkably, for the forward link the performanee differenee is higher than the one obtained 
in the return link. Note that, similar to the return link, the values of a are seleeted so that the 
feasibility of MSErl in (|3^ is hold. 

Einally, the dispersion index among users is analysed and represented in Eigure [H Eor this 
ease, the dispersion values are even higher for the referenee seenario and our approaeh leads to 
higher fairness between beams. 
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Fig. 6. Forward link throughput values. 



Fig. 7. Forward link availability. 
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Fig. 8. Forward link throughput index of dispersion. 



Fig. 9. Forward link throughput with respect to channel variations. 
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VII. Conclusion 

This paper proposes a design of non-ehannel adaptive beam generation proeess that inereases 
the system throughput eompared to the eonventional existing teehniques in both forward and 
return link of a multibeam satellite system. The design is based on an upper bound approximation 
of the worst ease SMSE, whieh results to be the same for both forward and return links, leading 
to a large reduetion of the payload eomplexity. The robust approximation relays on a first 
perturbation model whieh results tighter than eurrent robust designs. Moreover, the simulation 
results also have shown the potential advantage of the eonsidered design in order to inerease the 
total system throughput. As a eonsequenee, this new approaeh eould beeome a breakthrough in 
the design of the next satellite systems, whieh so far have designed the on-board beamforming 
only based on geographieal information. 


Appendix A 

The goal is to prove, the proposed optimal design of B in (1^ ean minimize the upper-bound 
of SMSErl in (iMl) . Eirst, by employing the eigenvalue deeomposition of Z in (1^ . problem 
(iMl) ean be rewritten as 


mm 


S.t. 


with the following definitions 


H \-l 


traee(Ii^ -f 


= Ik, 


Mfii = BtJ, 


and. 


Dh. a (s - 


0 


Kx(N-K) 


OrAT- 


{N-K)xK 0{N-K)x{N-K) 


(55) 


(56) 


(57) 


where S has only K non-zero eigenvalues, as has rank equal to K. Aetually, the problem 

(1551) ean be written as 


E 


1 

1 -I- 


mm 


(58) 
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S.t. 

where Ai(.) denotes the i-th largest eigenvalue of the respeetive matrix. Obviously, MDM^ is 
a hermitian matrix whose eigenvalues are always positive. Then, it follows that 

1 


9 {\) = 


i = 


1 + 

is eonvex funetion on XiCWLfn,!). By using the theorem 3.C.1 in ll24l . we have that 

1 


(59) 


K 


K 






( 60 ) 


i=l 


where A = (Ai(M/jLDijLMfi ),, and 0(.) is a schur-convex function 
operator. On other hand, the theorem B.l in lf24l proved that 

d ^ A, (61) 

where d(.) represents Kxl vector formed by the diagonal elements of the matrix 

i.e. d = ...,di^(MRLDRLMg^))^. Finally, combining of dM]) with the 

schur convexity of ^(.), we have that ^(d) < ^(A), i.e. 

1 1 


K 


E 


> 


E 


(62) 


^ i=i ^ + Xi(yLnL^RL^RL) 

Moreover, the equality in (l62l) is reached whenever is diagonal. To this end, it 

is clear that M has to be diagonal such that 

M/jl = [Ik 0kx{n-k)]- (63) 

Given (15^ . it implies that B has to be made of the K first rows of the matrix U^, that is 

B = L^.,k, (64) 

and concludes the proof. 

Appendix B 

The starting point of the derivation is the upper bound obtained when only considering the 
eigenvalues variation 

a (S - e^l) Uf, (65) 

where for this case we additionally consider the perturbation on the eigenvectors as 

(U, + AU,) (S - e^^l) (tJ, + AU,)^ . 


( 66 ) 
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In If20| it is presented that the perturbation on the eigenveetors take the form of 

AU, = U,R + u„u/azu,s;' 

where 


R = D ■ (Uf AZU,S + SUf AZ^U,) , 
and the g, /-th entry of D is 


for f ^ g and A/ for / = 1,..., A denote the eigenvalues of HH^. Considering that 


AZ < e^I, 


(67) 


( 68 ) 


(69) 


(70) 


the following inequality holds 

AU, < tj.R + e^UnUf (71) 

Additionally, we have that 

Uf AZU.S + EUf AZ^U, < e^Uf U.S + e^SUf U,. (72) 

The following lemma is required for obtaining the result 

Lemma 2 For any semidefinite positive matrices A, K, C, and K < C, it holds that 

A ■ K < A • C. (73) 

Proof. See Theorem 17 of fl251 . □ 

With this last result it is possible to write the following 

AU, < U,R + eHUnUn"'u.S;\ (74) 

where 

R = D ■ (e^Uf U,S + . (75) 
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